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Abstract 

In a supersymmetric extension of the standard model with local gauged baryon and lepton 
numbers (BLMSSM), there are new sources for lepton flavor violation, because the right-handed 
neutrinos and new gauginos are introduced. In BLMSSM, we study the charged lepton flavor 
violating processes Ij —>• Zj -|- 7 and Ij 3li in detail. The numerical results show that in some 
parameter space the branching ratios for charged lepton flavor violating processes can be large 
enough to be detected in the near future. 
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I. INTRODUCTION 


From the neutrino oscillation experiments [l|, it is convincing that neutrinos have tiny 
masses and mix with each other QQ. T.* .epto„ flavo. is „ot exact in 

the universe. Though the standard model(SM) has achieved great success with the detected 
lightest CP-even Higgs, SM should be extended. Because of the GIM mechanism, in SM 
the charged lepton flavor violating(CFLV) processes are very tiny, for example BrsMilj —t 
h + 7 ) The experiment upper bounds of the CLFV processes Ij li + j and 


Ij —)■ 3li are 


Br(/i —)■ ey) < 5.7 x 10 Br(/i —)■ 3e) < 1.0 x 10 

Br(r ey) < 3.3 x 10“®, Br(r /ly) < 4.4 x 10“®, 

Br(r 3e) < 2.7 x 10"®, Br(r 3/i) < 2.1 x 10"®. 


( 1 ) 


They are much larger than the corresponding SM theoretical predictions. To explore new 
physics beyond SM, study CLFV processes is an effective approach. Once physicists observe 
CLFV processes in future experiments, there must be new physics beyond SM. 

In a simple extension of SM, with a new additional Yukawa matrix for right-handed neutri- 

n 

nos, CLFV processes are induced at loop level with neutrino [6|]. They are suppressed strong 
by the tiny neutrino masses and impossible to be observed practically. One popular su¬ 
persymmetric extension of SM is the minimal suppersymmetric standard model(MSSM)[^. 
In R-parity conserved MSSM, the left-handed light neutrinos are still massless and can not 
explain the discovery of neutrino oscillations. Therefore, physicists extend MSSM to ac¬ 
count for the light neutrino masses and mixings. Adding low-scale right-handed neutrinos 
and approximate lepton number symmetries, z/^MSSM is obtained, where the authors study 
the CLFV processes 81]. In the supersymmetric standard model with right-handed neutrino 
supermultiplets, the authors investigate various LFV processes in detail ^ . In our previous 


work, we study neutrino masses and CLFV processes in /iz/SSMjinlj. 

For the beyond SM models, one can violate R parityjnj with the non-conservation of 
baryon number [B) or lepton number (L) 3, 13 j. A minimal sup ersy mmetric extension of 
the SM with local gauged B and L(B 


proposed in one of the references in 


/MSSM) is a favorite one 14j. BLMSSM was first 


14l |. In the work, this model is that we are adopting. 


The local gauged B is used to explain the matter-antimatter asymmetry in the universe. 
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Right-handed neutrinos are introduced in BLMSSM to account for the neutrino oscilla¬ 
tion experiments, which lead to three tiny neutrino masses through the seesaw mechanism. 
Then lepton number (L) is expected to be broken spontaneously around TeV scale. In 


BLMSSM, the lightes 


studied in the work 


CP-even Higgs mass and the decays —>■ 77 , —)■ ZZ{WW) are 

15| . Taking into account the Yukawa couplings between Higgs and exotic 


quarks, we study the neutron and lepton electric dipole moments(EDMs) in the CP-violating 
BLMSSM 3, 17|. mixing and t ^ c + ''y,t ^ c + g are also investigated in SM 

extension with local gauged baryon and lepton numbers [l^. 

In this work, we analyze these CLFV processes (/i ey, /r —)• 3e; r —)■ ey, r /ly, 
r —)■ 3e,r — )■ 3/i) in the frame work of BLMSSM. Compared with MSSM, there are new 
sources to enlarge these CLFV processes via loop contributions. The new CLFV scores are 
produced from 1 . the right-handed neutrinos mixing with left-handed neutrinos; 2 . the 
coupling of new neutralino(lepton neutrahno)-slepton-lepton. In some parameter space of 
BLMSSM, large corrections to the CLFV processes are obtained, and can easily exceed their 
experiment upper bounds. Therefore, to enhance these CLFV processes is possible, and they 
may be measured in the near future. 

After this introduction, we briefly summarize the main ingredients of the BLMSSM, and 
show the needed mass matrices and couplings in section 2. In section 3, the decay widths of 
these interested CLFV processes are analyzed. The input parameters and numerical analysis 
are shown in section 4 and we give our conclusion in section 5. 


II. BLMSSM 


BLMSSM is the supersymmetric extension of the SM with local gauged B and L, whose 
local gauge group is SU{3)c ® SU{2)i ® f/(l)y 0 V(l)s ® V( 1 )l 12 |. The exotic leptons 

L4 ~ (1, 2, -1/2, 0, L4), E 2 ~ (1, 1, 1, 0, -L4), N 2 ~ (1, 1, 0, 0, -L4), LI ~ 

(1, 2, 1/2, 0, —(3 -|- L 4 ))) F 5 ~ (1, 1, —1, 0, 3 -|- L 4 ) and ~ (1, 1, 0, 0, 3 -|- L 4 ) are 
introduced to cancel L anomaly. As well as, the exotic quarks Q4 ("NJ (3, 2, 1/6, R4, 0), 




(3, 1, -2/3, -H 4 , 0), Dl ~ (3, 1, 1/3, -R 4 , 0), Ql ~ (3, 2, -1/6, -(1 + 


B 4 ), 0), f /5 ~ (3, 1, 2/3, 1 -|- i? 4 , 0) and ZI 5 ~ (3, 1, —1/3, 1 -|- B 4 , 0) are introduced 
to cancel B anomaly. To break lepton number and baryon number spontaneously, the 
Higgs superfields <|)£,(1,1, 0, 0, — 2), (/)i(l, 1, 0, 0, 2) and <I>b(1, 1, 0,1, 0), </ib(1, 1, 0, — 1, 0) are 
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introduced respectively. Now, Higgs mechanism is the very massy foundtion stone for particle 
physics, and people are convinced of it, because of the detection of the lightest CP even Higgs 


at LHC 


19|. The Higgs helds and acquire nonzero vacuum expectation 


values (VEVs), then exotic leptons and exotic quarks obtain masses. In the BLMSSM, 
the superhelds X(l, 1, 0, 2/3 + H 4 , 0), X'(l, 1, 0, —(2/3 + H 4 ), 0) are introduced to make the 
heavy exotic quarks unstable. Furthermore X and X' mix together, where the lightest mass 
eigenstate can be a candidate for dark matter. 

The superpotential of BLMSSM is[f^ 

y^BLMSSM = y^MSSM + kV^ + kVi + kVx , 

kVs = XqQ4QI^b + ^uU^U^ifB + ^oDlD^ifB + ^^b^b'^b 

+Y^,Q^HuUl + Yd,Q4Hdbl + YuMHdU^ + Y^MHuD^ , 

Wl = Y,,UHdEl + Y.MH^Nl + Y^.LIH^E, + Y,,LIH,N, 

+YbLHuN^ + \n^N^N^0L + ^^L^L0L , 

kVx = XiQQlX + X2U%X' + X^b^b^X' + ^ixXX' . ( 2 ) 

where Wmssm is the superpotential of the MSSM. The soft breaking terms Csoft of the 
BLMSSM can be found in the works 14, 

^soft = - ml 02*02 

-m\fi^Ql - mlu;U, - ml DID, - mlV.U - mlNl*Nl 
-mlEl*E2 - mli^Ll - mlN;N, - mlElE, - 

- ml^ip*BipL - {rriBXBXB + miXLXi + h.c.) 
+{Au,Yu,Q4HuUI + Ad,Yd,Q4Hdb2 + Au^YuMHdU, + A^.Y^MHuD, 
AAbqXqQ^Q^^B + AbuXuU2U,(Pb + ABDXDb2b,(pB + OBfiB^B^B + h.C.j 
+ {A,,Y,,L4HdE2 + A,,Y,,L4H^N2 + A,,Y,,LIH^E, + A,,,Y,,LlHdN, 


YAxYj^LHuN^ + A}\fcX}sfcN'^N‘^(pi + BiUb^l^l + h.c. | 
+{HiAiQg^X + A2X2U^U,X' + A,X,b^b,X' + BxfixXX' + h.c.} . 

The SU{2)i doublets H^, should obtain nonzero VEVs Vu, Vd, 


( 3 ) 


/ 


H„ = 


H+ 


( 




Hd = 


X(^Vd\HlYiPt) 


HI 


( 4 ) 
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The SU{2)l singlets $ 5 , ips obtain nonzero VEVs vb, vb-, 

^B = ^{vB + ^% + iPl) , ^B = ^ivB + ^% + i^B) ■ ( 5 ) 

In the same way, the SU{2)l singlets obtain nonzero VEVs vl, vl, 

^L = ^{vL + ^l + iPl) , ^L = ^{vL + ^l + iPi) ■ ( 6 ) 

Therefore, the local gauge symmetry SU{2)l ® U{1)y <SU(1)b <S)U{1)l breaks down to the 
electromagnetic symmetry 17(l)e. 

= — sin + cos represent the charged Higgs, whose squared masses at tree 
level are . The charged Goldstone bosons and neutral Goldstone bosons 

ir± aO w ^ 

are denoted respectively 


= cos f3Hf + sin jSH^, = cos PP]I + sin (3P^ , 


T± 


T± 


Gl = cos 13gPi + sin = cos (3^Pi + sin P,, 


with tan/3 = tan/3g = Vg/Vg, tan/3^ = Vg/Vg. 


(7) 


= -sin/3P° + cos/3P°, = -sin/3,P° + cos/3,P„ , 


B B 


=-sin/3,P° + cos/3,P;', 


1,5 


( 8 ) 


are the physical neutral pseudoscalar fields, and their masses at tree level read as|l5| 

Pf^H _ 2 PbI^B _2 PL^^L 


m „ = 


m „ = 


m „ = 


(9) 


cos (3 sin (3 ^°b cos (3g sin (3g cos (3g sin (3g 

The lightest neutral GP-even Higgs is obtained from diagonalizing the mass squared 
matrix of neutral GP-even Higgs in the sector {Hi,Hi) 


( 




/ 


cos a sm a 


( 


— sm a cos a 


HI 

po 


tan 2a = — 5 ^ tan 2/3 . 
- I^AO 


( 10 ) 


and (p° mix together and the mass squared matrix is 


= 

EB 


^ m3 cos^/3„-1-sin^/3„, (m^ -1-) cos/3„ sin/3 

Zb \ Zg ^ _ 

{m3 + m3 ) cos/3„ sin/3„, m^ siv?l3„+m^ cos^/3 

^ ^B Zb ^ aO^ ^ 


\ 


( 11 ) 
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with VBt = \Jv'^g + "*^ 5 • '^zg — ds'^Bt denotes the mass of neutral U{l)g gauge boson Zg. 
Two mass eigenstates can be gotten 


/ 


B 


( 


( 


— sm ag cos 


$0 




( 12 ) 


by the mixing angle that is defined as 


m + m 
Zb ao 

tan 2a „ = — -r^tan2/3„. 

Zb 


(13) 


In the same way, we obtain the mass squared matrix for (p^) 


Kb = 


^ cos^ sin^ (5 ^, ) cos (3^ sin (3^ 

) cos (3g sin (3 ^, sin^ (3^ + rn?^^ cos^ (3^ j 


(14) 


with vl^ = 3- vK = 2mUf,. represents the mass of neutral U{l)g gauge boson Z^. 
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20 | analyze the mass matrices of exotic quarks, exotic squarks 


In BLMSSM, the authors 
and some exotic sleptons. 

With the introduced superhelds three neutrinos obtain tiny masses through the see¬ 
saw mechanism. After symmetry breaking, we obtain the mass matrix for neutrinos in the 
basis {u,N^)^ 

/ 


\ij 


(15) 


0 ^(n)^ 

Eq. ffTSj) can be diagonalized by the unitary matrix Uy. Then, one gets three light and 
three heavy neutrino mass eigenstates. 

The new gaugino \l and the superpartners of the SU{2^ singlets mix, which 

produce three lepton neutralinos in the base HXt . b’dij-. 'i/^,^)|l7|. 

^ 2Ml 2vLgL -‘^vlQl \ ( ih ^ 

‘^vlQl 0 -pL 
V -2vLgL -g-L 0 y 




L 

) 


+ h.c. 


(16) 


Using Z]\fg, one can diagonalize the mass matrix in Eq. fllbp to obtain three lepton neutralino 
masses. 
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In BLMSSM, the mass squared matrix of slepton is different from that in MSSM, because 
of the contributions from Eqs. (l2]l3il . The corrected mass squared matrix of slepton reads as 

/ 


{M\)ll {M\)lr 
\ {■^‘l)rr 

{MDll, {MDlr and {MDrr are shown here 


(17) 


(M 


lJll 


^IJ + gd'VL - Vl)Sij + miiSij + 


{M‘i)RR - 


gKvl-vj) 


{MI)lr = -+ -^iAi)ij, 


Sij - ain ~ "D^iJ + infiSrj + 


( 18 ) 


The unitary matrix is used to rotate slepton mass squared matrix to mass eigenstates. 

There are six sneutrinos, whose mass squared matrix is deduced from the superpotential 
and the soft breaking terms in Eqs. (l2]l3D . In the base = {9,N^), the concrete forms for 
the sneutrino mass squared matrix Aifi are shown here 

- d’u)^ij + ali^L ~ + (mi)ij, 

MKNTN'j) = -alivl - ii)S,j + f (PK)„ + 2«i(A]v.A;,.)„ 


— VuVl{Y^X] yc)ij H- -=(A]y)jj{Y^)jj. 


(19) 


The superhelds in BLMSSM lead to the corrections for the couplings existed in MSSM 
and some corrected couplings are deduced. We give out the couplings for W-lepton-neutrino 
and Z-neutrino-neutrino 

3 2 


^WLid — 
^Zvv 


\/2si4A 

e 


2swCiv I,J,K=liJ=l 

The charged Higgs-lepton-neutrino couplings are 


Z EiAdr'plrPLd^, 

I,J=1 i=l 

z, E EK.Y'ulUa’YPL 


Va 


( 20 ) 


3 2 


= E E cos l3UlhPL - FE* sin I3U‘^;uPr) E 

1^J=1 i=\ 


3 2 

-EE sin (XUIYPl + Yy* cos f5Ul\jPR) vi + h.c. 

I,J=1 i=l 


( 21 ) 
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The corrected chargino-lepton-sneutrino couplings read as 


3 2 


= -1; i; x-j[Y!^zT{zi„fPR 


I,J=1 i,j=l 



( 22 ) 


We also obtain the adapted Z-sneutrino-sneutrino couplings as follows 



(23) 


In BLMSSM, there are new couplings that are deduced from the interactions of gauge and 
matter multiplets — X^'ipiAj). After calculation, the lepton-slepton-lepton 

neutralino couplings are obtained 



(24) 


III. CHARGED LEPTON FLAVOR VIOLATION IN THE BLMSSM 

In this section, the CLFV processes Ij —)■ /* + 7 and Ij —)■ 3/* are studied in the BLMSSM. 
For convenience, the triangle, penguin and box diagrams are analyzed in the generic form, 
which can simplify the work. 

A. Rare decays Ij —)• k + l 

When the external leptons are all on shell, we can generally write the amplitudes for 

Ij —y "F 7 


M = ee>^Ui{p + q) Pz. + C^Pr) 

+ mlj^(y^,vq''{C2PL + C2PR)\uj{p) , 


(25) 


where p is the injecting lepton momentum, q is the photon momentum, and mi. is the mass 


of the j-th generation charged lepton. Ui{p) and Vi{p) are the wave functions for the external 
leptons. The relevant Feynman diagrams are shown in Fig{Tl The hnal Wilson coefficients 
Cf, C 2 , C 2 are obtained from the sum of these diagrams’ amplitudes. 






FIG. 1: The one loop diagrams for Ij ^ + 7 , with F representing Dirac(Majorana) particles. 


The contributions from the virtual neutral fermion diagrams Fig|T]^a) are denoted by 
a = 1, 2. We give out the deduced results in the following form, 


cHn) = E E 


F=x°S=L,L,H± 


FHp-HpFhixF.xs) , 


C^{n) = V V 

f=^xI,^s=Uh^ 

Cf(n) = Cf(n) a = 1,2, 


rrSFli rrSn^F 


l2{xF,Xs) - h{xF,Xs) 


(26) 


with X = m?/rri^ and m representing the mass for the corresponding particle. and 

Hl r are the corresponding couplings of the left(right)-hand parts in the Lagrangian. The 
one-loop functions Ii{xi,X 2 ),i = 1 ■■ - 3 are collected here 


Il{Xi,X2) = 


ll-|- 61 nx 2 15x2 + 18 x 2 In a ;2 , 6 x 2 -|- 18 x 2 In X 2 


967r2 L (x2 - xi) (x2 - XiY 

6 x^ In Xi — 6 X 2 In X 2 - 


+ 


(X2 -Xi) 


3 ’ 


h{xi,X2) = 


(X2 - Xi) 

3 -|- 21 nx 2 2 x 2 -|- 4 x 2 In X 2 2x^lnxi 2 x 2 lnx 2 ■ 

327r2 L (X2 - Xi) (X 2 - Xif (X2 - Xi)^ (X2 - Xi)^- 


N 1 rl-|-lnx 2 Xi Inxi — X 2 lnx 2 i 
h[xi,X2) = 7 -r + 


(27) 

(28) 

167r2 L(a;2 - Xi) ' (^X 2 -Xif 

CY^'^{c),a = 1,2 stand for the coefficients from the virtual charged fermion diagrams 
FiglH^b), and they are shown here 

Ci'(c) = ^'^\h{xF,xs) - 2I4 {xf,xs) - Ii{xf,xs) 

bm 


F=x^ S=i^ ‘‘W 

ci{c)= Y. E 




C^ic) = C(c 


f=y± s=i> ^ij^w 

L-tr^Ri a = 1 , 2 


h{xF, xs) - h{xF, xs) - Ii{xf, xs) 


=X-‘ 

-iLi 


with 


/4(a;i,X2) = 


l-|-lnxi Xi Inxi — X 2 lnx 2 


IGtt^ 1 (X2 — Xi) 


(X 2 - X,Y 


(30) 


(31) 


9 









































Because three light neutrinos and three heavy neutrinos mix together, the virtual W 
diagrams FigJT]^c) have corrections to the CLFV process Ij —)■ lij. The corresponding coef- 
hcients are denoted by C^’^(hF)(Q; = 1 , 2 ) 


WFli 




F=v 


h{xF,Xw) + h{xF,Xw) 


C^m = Y. + ^)\2I^{xf^xw) - \h{xF.xw) 

C'^(hF) = 0 , a = 1 , 2 . 

The total coefficients are the sum of Eqs. (l26lhl30lhl32|) 

CL,r ^ + C^^ic) + Ct^iW), ^ = 1,2. 

With the Eq. (l25D . the decay width for Ij k + j can be expressed as 




a 


+ 


a 


R 


(32) 


(33) 


(34) 


B. Rare decays Ij 3 h 


The CLFV processes Ij —)■ 3li are very interesting. Both penguin-type diagrams and 
box-type diagrams have contributions to the effective Lagrangian. With Eq.(|25]), one can 
obtain the y-penguin contributions in the following form. 


T^-p = Ui{pi) g^7^(C'fPL + C^Pr) + mi.iap^q'^iC^PL + C^Pr) Uj{p) 


Ri 


X -^Ui{p2)^^Vi{p3) - (Pl -H- P 2 ) • 


(35) 



kipi) kip?) 


FIG. 2: The penguin-type diagrams for CLFV process L 


3li. 


The contributions from Z-penguin diagrams are depicted by the Fig 12], and deduced in 
the same way as 7 -penguin diagrams. 


Tz-p = 


m‘ 


-u,{p{)^p{NlPl + NRPR)u^{p)u,{p2)r{Hl^^'^PL 
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( 36 ) 


+ - (Pl ^ P2) , 

Nl,r = Nl,r{S) + Nl,r{W) . 

The concrete forms of the effective couplings Nl{S), Nr^S) read as 


2771 ^, 771^2 ^SF 2 k TTZF 1 F 2 uSnjF,_^ ^ ^ \ 

2 ^R ^l\XSi XP 21 Xp^) 


m 


w 


^ i: i: 

^=X°,X^,!'S=L,l>,B± 

+ E i: [//?"■ la.la 

F=xlS=L 

Nr{S) = iVt(S)U„H. 


(37) 


The functions Gi{xi,X 2 ,X 3 ) and ^ 2 ( 2 ^ 1 , 2 ^ 2 ,a^s) are 

1 r xilnxi 


Gi(a;i,a;2,a;3) = 


+ 


X2 In X2 


16772 L {xi - X2){xi - X3) (X2 - Xi)(x2 - X3) 

X 3 In X 3 


G2{Xi,X2,X3) = 


{X3 -Xi){x3-X2)^' 


(A + 1 + Inx^) + 


(38) 


xf Inxi 


16772 L 

+ 


Xn In X 2 


+ 


{xi - X2)(xi - X3) 

xl In 0:3 


((C2 - a;i)(a;2 - 0 : 3 ) ( 0:3 - a(i)(a;3 - 2 : 2 )^ 


(39) 


G 2 {xi,X 2 ,X 3 ) has inhnite term, and to obtain hnite results we use MS subtraction and DR 
scheme. 

We deduce the effective couplings Np p(W) in detail and keep the small terms. 

Nl{W) = —^ G3{xp,xw) +‘2{xi + Xj)[Ii{xp,xw) — h{xp,xw)] 

esw p=y 

_j_ \ ' j^WF 2 li jjW*ljFi j^z*F1F2 ^ ^ 

' 2 / ^ L L L 


Fi^F2=D 


( - ~ ^‘2(^W,XPJ^,XP2) 


+Xj[-G4{xw: Xp^^Xp^) + G3{xw, Xp^^^Xp^)]^, 


Cw 


Nr{W) = \ 2 ,/^[h{xp,xw) - l2{xp,xw)] 

P = iy 

1 


E 

Fi,F2=D 


WF2li jjW*ljFi jjZ*FiF2 


y/XiXj (2G1 (xw , Xp^ , XP2 


—-Gi{xwiXp^^xp2) — 2G^{xw-,xp^^xp2)^. 


(40) 


The concrete expressions for the functions G 3 (a;i, 0 : 2 ), G 4 ((ri, 0 : 2 , 0 : 3 ) and G^{xi,X 2 iX 3 ) are 
collected here 

-1 

l(A + lna;^ + lj + — 

16772^ (2:2-a;i)" 


G3(2;i,2;2) = 


/ xl\nx2-xl\nxi X2 + 2x2lnx2 1\ 

( A + Inx^ + 1 +- j - 7 ^-+-77 ), 

^ I.Tn — Til"! 2 ;^ _ 2 ;^ 22 
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G4.{xi,X2,X3) 


G5{xi,X2,X3) 


1 / 2 a;f [ 3 a;i(a;i — X2 — X3) + X2 + X2X3 + x^] Inxi 

327r2 ^ {xi — X2 Y{xi — X 3 Y 

2 ( 3 a;f — ?iXiX2 + X2)x2 In X2 2 ( 3 a;^ — 3a;ia;3 + x'^x^ In 0:3 

{xi - X2 Y{x 2 - X3) (xi - X3)^(X2 - X3) 

Xi [bxj - 3 xi{x2 + X3) + X2X3\ \ 

(xi-a;2)2(a;i-0:3)2 /’ 

1 /xf (2a:i — 0:2 — 0:3) Inxi a:2(a:2 — 2a:i) lna:2 

(^xi - X2y{xi - X3y {xi - X2y{x2 - X3) 

__^ 3:3(23:1 - X3) \nx3 X 

(3:1 - 3:2)(3:1 - 3:3) (3:1 - a:3)2(3:2 - 3:3)^' 


(41) 
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hip) 

F k(pi) 

hip) 

F kipi) 

hip) 

kipi) 



1 
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s', 

is 

si 
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si 

tfe); 

; kipi) 

t(P 2 ) 

kiPi) 

kiP2) \y^ ^ 

\ kipa) 


F 


F 




(a) 

hip) F 

kipl) 

(b) 

hip) 

(c) 

kip\) 





* >\F 




hw 

wi 

0 i> 




(i(P 2 )> 

yiipi) 

kiP 2 )Sx 

^\^kiPi) 



F 

(e) 


(f) 


FIG. 3 : The box-type diagrams for CLFV processes Ij 3 li with F representing Dirac (Major ana) 
particles. 

The box-type diagrams drawn in Figj3]can be written as 

Tbox = {B^e‘^Ui{pi)j^PLUj{p)ui{p2)'j^PLVi{p3) + (T O R)^ 

+ [p^e^ Ui{pi)'y^PLUj{p)ui{p2)j^PRVi{p3) - (pi O P2) + {L -H- /?)} 

-h [p^e^ Ui{pi)PLUj{p)ui{p2)PLVi{p3) - {pi ^ P2) + (T i?)} 

-h [P^e^ Ui{pi)af,uPLUj{p)ui{p2)a'^‘'PLVi{p3) - (pi ^ P2) + (T f-)- i?)| 

-h [P^e^ Ui{pi)PLUj{p)ui{p2)PRVi{p3) - (pi O P 2 ) + (T o R)Y (42) 


From the box-type diagrams, we obtain the virtnal chargino contribntions to the effective 
conplings with /3 = 1... 5 

1 


bF)= E E 

Fi,F2=x^ Si,S2=v 


Bi(c) = E E 

Fi,F2=x^ Si,52=!> 
Tflp-^Tflp^ 


(rr rr rr rr \ TjF2Flli TjS\ljF\ prSlF2li prS2liF2 

^F2i ^Sii ^S2 )Pr Pr Pr Pr ) 


n (rr rf. rr ^ TJS 2 FRi ttSiI jFi RiF2li TTS2liF 2 

^&\FFi) ^F 2 ) ^S\) ^S 2 )PR Pr Pr Pr 




2 e2m^ 


(rr> ^ ^ 'y ) J-T^‘ 2 ^lh TTSiF 2 li JjS 2 liF 2 

XF 21 ^ 82 )^R 


Bt(c) = E E 


'Ulp-^Tflp^ 

Fi,F2=x± Si,S2=S ^ '^W 
?R 


(rr ^ -T* -T* ^ Jj 82 F\li JPrSlF 2 li TTS 2 liF 2 

(j'jyXFi"! ^F2'! ^Si"! ^S2)L ’ 


Bi(c) = BPc) = 0 , B^(c) = B^c) /? = 1... 5 . 


( 43 ) 


with 


GQ{Xi,X 2 ,X 3 ,Xi) = 


+ 


x\ Inxi 


+ 


X 2 In X 2 


167 r 2 L {xi - X2)ixi - X3) (xi - X4) (X2 - Xi)(x2 - X3) (X2 - X4) 


Xn In X3 


+ 


xl In X 4 


{x3 - (ri)(a;3 - X2){x3 - X 4 ) {X 4 - a;i)(a;4 - X2)ix4 - Xg)^ ’ 
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Gy{Xi,X2,X3,X4) = 


X\ Inxi 


+ 


X 2 In X 2 


167r2 L (xi - X 2 )ixi - Xs) (xi - X4) {X 2 - Xi){x 2 - X3)(X2 - X4) 

X 3 In X 3 X 4 In X 4 


+ 


( 44 ) 


(Xs - Xi)(x3 - X2)(X3 - X4) {X4 - Xi){x4 - X2){X4 - ' 

For the box-type diagrams, the neutralino-slepton, neutrino-charged Higgs and lepton 
neutralino-slepton contributions to the effective couplings are gotten, 


Bi(n) = y 
2 = 

1 


E 


Tflp^Tflp^ ^ I ^ ^ ^ ™ ^ ZTS2F\li TjS2F2li o-‘S'pi^2 


F,,F,=x^,xlx SpS 2 =l:l,H± 


'■ n ( rf rr. ry ry \ TT ^2 B lli TT1''j ^ 1 pf i>2 B 2li TTI 

-Lj^yXp^^ XP 21 XXS 2 )^L 


2 e^rn^ 


Gq i^xp^ , X P 2 5 Xp^ , X ^2 ) 


TTS 2 Flli TjSlljFl jjSlF 2 li TT^ 2 ^iF 2 

R L R L 


J^S 2 F\li jj^S*ljF\ j^g 2 F 2 li 
L R R L 


Bi(n) = 




Fi,F2=x°,xl,i' Si,S2=L,L,H± 


Tnp-^mp2 . \TTS2F1li pr^lhFl prSlF2li TTF2kF2 

n 2 4 ^ 7 \^Fii ^F2 i ^Sii ^82)^R ^R 


+ 


4e^mM 


'G^ {^Xp^ , Xp 2 1 Xg^ , Xg 2 ) 


w 


TJ-S2Flli TjSlljFl jjSlF2li TT^2^iF2 

R L L R 


prS2F\li TjS2F2li pr^t^i^‘2 

^R ^R 


B^(n) = 




Fi,F2=x°,xI,>^ Si,S2 =L,L,H± ^ 

^ j^S 2 FiIi j^S^ljFi jyS 2 F 2 li j:jS*liF 2 

2 ^ L L L 


mp^mF 2 ( ^ 

2 4 G'jyXp^^Xp2-,Xg^^Xp2) 


jjS 2 Fili IjFi jjSiF 2 li jjS 2 kF 2 


B^(n) = 


Fi,F 2 =x°,x 1 ,‘'Si,S 2 =L,L,H± ^ 


B-(n) = 0, B^(n) = B-{n] 


L^R, /3 = 1...5. (45) 


We also deduce the box-type contributions from virtual W-neutrino 


bHw)= E 

Fi,F2 

mp^mp 2 d 


1 


P,,P,=u 


d 


-2 


dxw 


dxw 

Gi{xw,xp^,xp 2 )H^ ^ 




WljFi ^w*F2Tij^WliF2j^W*Fji 
L L 


f>L(xxx^ sr ^ rnp,mp 2 d 

3 (^) z _^ 2 e^m*^ dxu / 

p^ p^—y ^ C nl^ UXm 


g4{x^, xp„xp 2 )HY^^^xr^^^^Hr^^Hr^^\ 
1 


B^{W) = 0, B^{W) = -B^{W), B^{W) = --B^{W), 


Hf (W) = B^{W) = 0, H"(W) = BpW), t = 3,4, 5. 




?R 


(46) 


With Eqs.( 


yL,R 


are expressed as 


bL,R ^ ^ ^ i?^’"(W), (/3 = 1... 5) . 


jL,_R 


L,R/ 


(47) 
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The decay widths for L —)■ Sk can be computed from the front amplitudes, 


r(/j —>■ 3/j) — 


-m'; 


5127r3 


{( 


ci 


+ 


3 2mz, 9 ' 


+ ( 


C 


1 

1 , _r 2 




^2 
yL/^R* 


) - 2(C'fC'f + + H.c.) + -( 

6 


+ 


R 


Bi 


B!, 


B. 


R 


R 


1 


Bi 


+ 


5 




^) + -(53^ + Bi ) + 6(5f + 5. 

) - ^{B^B^* + + C^B^* + + H.c 
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+ \{CiB^* + + H.c.) - ^NlrB^* + NRLBi* + H.c.) 

3 0 

- ]^{B^Bi* + + H.c.) + + C'fHf* + C^B!^* 


yR d-R* 

1 r 


Cfaf + H.c.) - -(C£Bt"‘ + C-jBr + C^B^ + C^B“ + H.c.) 


3L 


2(\Nll\‘ + lAfREl") + (liViRr + liVsiR + [BiNR + 


+ Bi_Nl^ + + H.c.) + 2(C^NR + + H.c.) 

+ (CfAfI„ + CfiVJi + H.c.) - MC«NR + C^-AfJs + H.c.) 
-2(C^/VJi + CfAf2„ + H.c.)]}, 


(48) 


with 


Nll = 


Nrr 

= Nrr 

2 ’ 
m^z 

Nlr = 

9 5 

Nrl 

= Nrr 


mt 


With the fomula the branching ratios of Ij Sk are obtained. 


(49) 


IV. NUMERICAL RESULTS 


In this section we discuss the numerical results, and consider the experiment constraints 
from the lightest neutral CP-even Higgs mass ~ 125.7 GeV lOj] and the neutrino 
experiment data. In this model, the neutron EDM, lepton EDM and muon MDM are 
studied in our previous works, and their constraints are also taken into account. In this 
work, we use the parameters |2lt] L 4 = |,A 7 vc = 1- The Yukawa couplings of neutrinos 
{YyY'i (/, J = 1, 2, 3) are at the order of 10“® ~ 10“®, whose effects to the CLFV processes 
are tiny and can be ignored. 
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To simplify the numerical discussion, we use the following relations 


iAi)ii = AL, {ANc)ii = {AN)ii = AN, 

= s^, for i = 1,2, 3. (50) 

If we do not specially declare, the non-diagonal elements of the used parameters should be 
zero. 

A. Ij —y h T 7 

Charged lepton flavor violation is related to the new physics, and the branching ratio of 
the process /i ^ e-l-y is strict. Its experiment upper bound is 5.7 x 10“^^ at 90% confidence 
level. At this subsection, the supposed parameters are AN = —500GeV, = 3TeV. 

1 . /i —)■ e - 1 - 7 

With the parameters Mgn = lTeV,/iL = neV,gL = 1/6, VLt = 3TeV, tau/dx, = 2, we 
numerically study the CLFV process /r —>■ e - 1 - 7 . The mass matrix of neutralino includes mi 
and m 2 , m 2 is also related with the mass matrix of chargino. Therefore, the two parameters 
mi and m 2 can effect the contributions(neutahno-slepton, chargino-sneutrino) for /r —)■ e-l -7 
to some extent. Supposing g,H = 480GeV, tan/3 = 11, Sm = 3400GeV, AL = —2TeV, = 
300GeV, we scan the parameters of mi versus m 2 in FigJU It implies that in this condition 
mi should be in the region (430 ~ 630)GeV and the effects of m 2 are small. 

1400 

1200 

1000 
d 800 

S 

600 

400 

200 

500 1000 1500 2000 2500 3000 3500 4000 

m2/GeV 




FIG. 4: For p —?■ e -|- 7 , the allowed parameters in the plane of mi versus m 2 with gn = 
480GeV,tan/3 = 11, Sm = 3400GeV,AL = — 2TeV, A)^ = SOOGeV. 

The slepton mass squared matrix has A/ and AL as non-diagonal elements which affect 
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the results through slepton-neutralino and slepton-lepton neutralino contributions. With 
= 480GeV, tau/S = 12, Sm = 3300GeV, mi = SOOGeV, m 2 = ITeV, in FigJSjWj;^ versus 
AL are scanned. FigjS] shows that the effects from AL are smaller than the effects from 
The allowed region of A'j^ is about (—2 ~ 2)TeV. 

4000 

2000 

> 

<u 

% 0 

-2000 

-4000 

-4000 -2000 0 2000 4000 

AL/GeV 

FIG. 5: For p —?■ e + 7, the allowed parameters in the plane of AL versus A'j^ with fin = 
480GeV,tan/3 = 12, Sm = 3300GeV,mi = 500GeV,m2 = ITeV. 

tan/5 is related to the mass matrices of chargino, neutralino, slepton and sneutrino, espe¬ 
cially to the non-diagonal elements of these matrices. In BLMSSM, almost all contributions 
to GLFV processes are influenced by tan /5. It is a sensitive parameter and affects the numer¬ 
ical results forcefully. m| and m^ are introduced in the soft breaking terms. Both slepton 
and sneutrino mass squared matrices include m| and m^ , which should give considerable 
effects to GLFV processes. Supposing hh = 470GeV, mi = SOOGeV, m 2 = ITeV, AL = 
—2TeV, 4^ = 300GeV, we plot the results with the allowed tan/5 versus Sm in FigJOl As we 
expected, they both are sensitive parameters. Because the upper bound of i?r(/i —)■ e - 1 - 7 ) 
is small, it is easy to exceed the bound in BLMSSM with the new contributions. 

2 . r ^ e -|- 7 

The experiment upper bound of i?r(r —)■ e - 1 - 7 ) is 3.3 x 10“® which is almost five-order 
larger than that of Br{fi e -|- 7 ). Here we use the parameters m 2 = ITeV, tan/5 l = 
2,Msn = 2TeV, 4L = — 2TeV, 4^ = 300GeV,= ITeV, (/i = 4 and vl^ = 3TeV. As 
discussed in the previous part, Sm can affect the contributions strong through slepton masses. 
Both slepton-neutralino and slepton-lepton neutralino give one loop corrections to the GLFV 
processes. Using the parameters mi = 600GeV,/in = 700GeV and tan/5 = 10(15,25) 
in Fig. ([7]) we plot the results varying with Sm by the dashed line, dotted line and solid 
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FIG. 6 : For fx —)■ e + 7 , the allowed parameters in the plane of tan 13 versus Sm with fXH = 
470GeV,mi = 500GeV,m2 = lTeV,^L = -2TeV,^i = 300GeV. 

line. These three lines are all decreasing fnnctions of the enlarging Sm- In the Sm region 
(1000 ~ 1500)GeV, the dashed line varies from 1.0 x 10“® to 1.0 x 10“®; the solid line varies 
from 1.0 X 10“^ to 1.0 x 10“®. As Sm > 2500GeV, the three lines are all much smaller than 
the upper bound. Gorresponding to same Sm during the region (1 ~ 2)TeV, the solid line 
results are about 10 times as the dashed line results, and the dotted line results are about 3 
times as the dashed line results. It implies that both Sm and tan (3 are sensitive parameters 
to the numerical results. Larger tan (3 leads to larger results obviously. 



FIG. 7: With tan/3 = 10(15,25), the results of Br{T ^ e + 7 ) versus Sm are plotted by the dashed 
line, dotted line and solid line respectively. 

fiH is included in the mass matrices of chargino and neutralino, which should produce 
considerable influence on the numerical results. With small tan/3 and large Sm, the results 
for r —)■ e + 7 are much smaller than the experiment upper bound. To embody effects from 
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f^H and mi, we use large tan/9 = 25 and small Sm = SOOGeV. The allowed numerical results 
are plotted by the dots in the FigEl 
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QJ 

S. 0 
1 ; 
a, 

-2000 

-4000 

-4000 -2000 0 2000 4000 

mi/GeV 

FIG. 8 : For r e -|- 7 , the allowed parameters in the plane of versus mi with Sm = 
SOOGeV, tan/3 = 25. 

The non-diagonal elements of (m|) and represent the lepton flavor mixing and 

lead to strong mixing for sleptons (sneutrinos). To simplify the discussion, the assumption 
{jnY)ij = = MLf for i ^ j and = 1,2,3 is used. We also consider the non¬ 

diagonal elements of Ai with the supposition (4;)^- = Af for i ^ j and = 1,2,3. Using 
the parameters hh = 480GeV, mi = SOOGeV, = ISOOGeV, tan/3 = 15, in the plane of 
MLf versus Af the parameter space is scanned, and the allowed results are shown in FigJO] 
The effects from MLf are stronger than the effects from Af. Here, the allowed region for 
MLf is about (—4 x 10® ~ 8 x 10®) GeV^. 
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FIG. 9: For r e -|- 7 , the allowed parameters in the plane of MLf versus Af with = 
480GeV,mi = 800GeV,Sm = 1500GeV, tan/3 = 15. 
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3. r —)■ /i + 7 

Similar as r ^ e + 7 , the branching ratio of r ^ /i + 7 is also large, whose experiment 
npper bonnd is 4.4 x 10“®. For the decay r ^ /i + 7 , the parameters m 2 = neV,fiH = 
SOOGeV, mi = SOOGeV, = 2TeV, AL = —2TeV, Sm = ITeV, A'j^ = 300GeV are nsed. 
The gangino mass mi is related to the nentralino-slepton contribntions to GLFV processes. 
With tan/3 = 15, hh = 500(1500, 3000)GeV, in Fig. ffTOl) the resnlts of Br{T —)■ /x + 7 ) 
are stndied versus mi by the dashed line, dotted line and solid line. As |mi| is around 
OOOGeV, these three lines arrive at their big values. The results are decreasing functions of 
the increasing |mi|, when |mi| is larger than 800 GeV. The biggest value of the dashed line 
can reach 3.2 x 10“^. The solid line varies from 1.0 x 10“^° to 5.1 x 10“®. The dashed line 
is the highest line and the solid line is the lowest one. The dotted line is the middle line 
and at the order of 10“®. The three lines show the GLFV processes are suppressed by heavy 
virtual particles at several TeV order. 



FIG. 10: With fin = 500(1500, 3000)GeV, the results of Br{T —p + 7 ) versus mi are plotted by 
the dashed line, dotted line and solid line respectively. 

Gompared with MSSM, tan /3l and are new parameters that have relation with mass 
matrices of slepton, sneutrino and lepton neutralino. Therefore, the effects to GLFV process 
r —)■ p + 7 from tau/di, and are of interest. Based on the supposition gi = = 

ITeV, tau/d = 18, we scan the parameter space of tau/d^ versus VLt in FigJTTJ The value of 
tan(3L can vary from (0 ~ 40), whose effects are small. As tan/^^ > 2 . 0 , vit should be no 
more than 3600 GeV. 

ql is not only the coupling constant for lepton neutralino-slepton-lepton but also con¬ 
stitute the mass matrix of lepton neutralino. Gonsiderable influence to r ^ p - 1 - 7 from 
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FIG. 11: For r —>•// + 7 , the allowed parameters in the plane of tan/^L versus vit with gi = 
= lTeV,tan/3 = 18. 

Ql is hopeful. The new gaugino mass fii is the non-diagonal element of the lepton neu- 
tralino mass matrix. To see how gi and /i^ affect the numerical results for r ^ /i -|- 7 , with 
ta.n(3L = = 3TeV, tan/3 = 18 we give out the allowed dots in the plane of ql versus 

fiL- FigJT2] implies that when gi is near 0.5, the results are larger than the experiment upper 
bound. The effects from /Xi is very weak, and can be neglected. 
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FIG. 12: For r ^ /i + 7 , the allowed parameters in the plane of gi versus gi with tan/?/, = 
2,r^Lt = 3TeV,tan^ = 18. 


B. Ij ^ 

In this subsection, we numerically study the CLFV processes Ij Sk with the supposed 
parameters mi = 3TeV,= lTeV,x;Lt = 3TeV. These processes have close relations to 
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4 + 7 - 

1. /i —3e 

The most strict branching ratio of CLFV processes Ij —)■ 3/j is Br{^ —)■ 3e), whose 
experiment npper bonnd is 1.0 x 10“^^. This experiment constraint is the first one to be 
considered for Ij —)■ 3/j. To stndy the process /i —)■ 3e, the nsed parameters are = 
ITeV, Sm = 3300GeV, AiV = — SOOGeV, AL = — 2TeV, = 300GeV, mi = 500GeV, (7^ = 
1/6, tan= 2. From the discnssion in the front section, tan/9 and non-diagonal elements 
of (m|) and are sensitive parameters to the GLFV processes. With = SOOGeV and 
7712 = ISOOGeV, the parameters tan/9 versus MLf are scanned in the Fig{T3l The plotted 
dots represent the allowed results which embody the influences from tan/9 and MLf. Here, 
the value of tan (3 should be no more than 15. 
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-20000 0 20000 40000 60000 
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FIG. 13; For p — 3e, the allowed parameters in the plane of tan/3 versus MLf with = 
500GeV,m2 = ISOOGeV. 

Here we consider the non-diagonal elements of (m|) and and suppose MLf = 

lO^GeV^ and tan/9 = 10. After the numerical calculation, the allowed parameters in the 
plane of m 2 versus fin are shown in the FiglTTl When hh is near 500 GeV, m 2 can vary from 
-3 TeV to 3 TeV. As > 600 GeV, the allowed scope of m 2 shrinks with the enlarging ^h- 

2. r ^ 3e 

The experiment upper bound of the GLFV process i3r(r —)■ 3e) is (2.7 x 10“^), and 
it is about four order larger than that of /i — )• 3e. For the study of r ^ 3e, tan {3 = 
10 ,Ms„ = 2TeV, fiH = 3TeV, m 2 = 1500GeV, mi = 2500GeV are supposed here. To show 
the importance of the non-MSSM contributions from lepton neutralino-selepton, we discuss 
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FIG. 14; For yu 3e, the allowed parameters in the plane of m 2 versus fin with MLj = 
10^GeV2,taii/3 = 10. 

the effects from ql and tan/3i with Sm = 3500GeV, AiV = —SOOGeV, AL = = 

300GeV. Fig {15] implies in the ql region (0 ~ 2), the allowed scope of tan/^^ is from 0 to 50. 
When ql is larger than 2.2, the region of tan/3i tnrns very small which is just from 0 to 2. 
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FIG. 15: For r ^ 3e, the allowed parameters in the plane of tan/^^ versus gi with Sm = 
3500GeV, AN = -500GeV, AL = -2TeV, A'^ = 300GeV. 

(m|) and (m^) are sensitive parameters relating to lepton mixing between different gen¬ 
erations. That is to say, their diagonal and non-diagonal elements are all important factors 
for GLFV processes. As AN = 2TeV, AL = 2TeV, = 500GeV, 5 :^ = 0.1, tan= 1, 
the allowed scope of Sm versus MLf is plotted in the Fig{T 6 l MLf should be no less than 
—2.0 X lO^GeV^ and the allowed smallest values of MLf turn large with the enlarging Sm- 
3. T ^ 3fi 
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FIG. 16: For r —>■ 3e, the allowed parameters in the plane of Sm versus MLf with AN = 
2TeV, AL = 2Te\,Ai^ = SOOGeV,^^ = 0.1, tan/3£, = 1. 

Similarly, we calculate the CLFV process r —)• 3p, whose experiment upper bound is 
Br{T —)■ 3/r) < 2.1 x 10“®. To obtain the numerical results for r —)■ 3/i, we use Mgn = 
lTeV,p// = SOOGeV, m 2 = ISOOGeV, mi = 2500GeV, AiV = 3TeV, AL = 3TeV, = 
3 TeV,5ii, = 0.4. In the studied processes, there are two angles tau/S and tau/dx, relating to 
the contributions. In the plane of tan/? versus tan/5x, with Sm = ISOOGeV we show the 
allowed results denoted by the dots in Fig{T71 The suitable value of tan/3 is in the region 
(0 ~ 10). Gompared with the effects from tan/3, those effects from tan/3x are tiny. 

30 
25 
20 
t 15 
10 
5 
0 




FIG. 17: For r —>■ 3/x, the allowed parameters in the plane of tan/3 versus tan/3x with Sm = 
ISOOGeV. 

The mass squared matrix of sneutrino include An, An^ and which naturally influ¬ 
ence the contributions from sneutrinos and charginos. We take into account the non-diagonal 
elements of these parameters and suppose {AN)ij = {AN<^)ij = ANf, = MNf for 
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i ^ j and i,j = 1,2,3. With tan/? = 10 , tan= 10, Sm = 2TeV, in FigJTHlwe plot the 
allowed results in the plane of ANf versus MNf. To obtain the allowed results, ANf is no 
more than zero. The effects from MNf are tiny and ignorable. 
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FIG. 18: For r —>■ 3p, the allowed parameters in the plane of ANf versus MNf with tan/3 = 
10, tan Pl = 10, Sm = 2TeV. 


V. DISCUSSION AND CONCLUSION 

In SM the theoretical predictions for CLFV processes Ij ^ + 7 and Ij —)■ 3li are much 

smaller than their experiment upper bounds. If large branching ratios of CLFV processes 
are detected, there must be new physics beyond SM. In BLMSSM, there are new parameters 
and new contributions to CLFV processes. For example, beside three light neutrinos there 
are three heavy neutrinos and six sneutrinos. Furthermore, new gauginos and new higgsinos 
mix leading to three lepton neutralinos, that can give new type contributions through the 
coupling of lepton neutralino-slepton-lepton. 

The branching ratio experiment upper bounds of /i e + 7 and /i —)■ 3e are strict 
and rigorously conhne the parameter space. For the both processes, it is very easy to 
exceed their experiment upper bounds. The experiment upper bounds of the processes 
r—)'e + 7 ,r—)-/i-|- 7 ,r—)-3e and r —)■ 3/2 are all at the order of 10“® and much larger 
than those of /i —)■ e + 7 and fi —)■ 3e. In our used parameter space of BLMSSM, the 
branching ratios of these six CLFV processes can be large enough to achieve the bounds 
and even surpass them. From the numerical results, one hnds the important parameters are 
tan (3, Sm, MLf, mi, m 2 , and g^, where tan/?, Sm, MLf are very sensitive parameters. We 
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hope in the near future large branching ratios of CLFV processes can be detected by the 
experiments. 
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